Interactions and Interference in Quantum Dots: 
Kinks in Coulomb Blockade Peak Positions 
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We investigate the spin of the ground state of a geometrically confined many-electron system. 
For atoms, shell structure simplifies this problem- the spin is prescribed by the well-known Hund's 
rule. In contrast, quantum dots provide a controllable setting for studying the interplay of quantum 
interference and electron-electron interactions in general cases. In a generic confining potential, the 
shell-structure argument suggests a singlet ground state for an even number of electrons. The 
interaction among the electrons produces, however, accidental occurrences of spin-triplet ground 
states, even for weak interaction, a limit which we analyze explicitly. Variaton of an external 
parameter causes sudden switching between these states and hence a kink in the conductance. 
Experimental study of these kinks would yield the exchange energy for the "chaotic electron gas" . 



The evolution of the properties of a system as a con- 
tinuous change is made to it is a ubiquitous topic in 
quantum physics. The classic example is the evolution 
of energy levels as the strength of a perturbation is var- 
ied . Typically, neighboring energy levels do not cross 
each other, but rather come close and then repel in an 
"avoided crossing". However, if there is an exact sym- 
metry, neighboring levels can have different symmetries 
uncoupled by the perturbation, and in this special case 
they can cross. 

A new and powerful way of studying parametric evo- 
lution in many-body quantum systems is through the 
Coulomb blockade peaks that occur in mesoscopic quan- 
tum dots The electrostatic energy of an ad- 
ditional electron on a quantum dot- an island of con- 
fined charge with quantized states- blocks the flow of 
current through the dot- the Coulomb blockade [^,0. 
Current can flow only if two different charge states of 
the dot are tuned to have the same energy; this pro- 
duces a peak in the conductance through the dot. The 
position of the Coulomb blockade peak depends on the 
difference in ground state energy upon adding an elec- 
tron, Eg^,.{N+l)-Eg^{N). Thus, the evolution of the 
peak position reflects the evolution of these many-body 
energy levels as a parameter, such as magnetic field or 
gate voltage, is varied. 

Since quantum dots are generally irregular in shape, 
the orbital levels have no symmetry and so avoid cross- 
ing. However, the spin degrees of freedom are often de- 
coupled from the rest of the system, so states of differ- 
ent total spin can cross. Such a crossing will cause an 
abrupt change in the evolution of the Coulomb blockade 
peak position- a kink- as the spin of the ground state 
changes, even though there is no crossing in the single- 
particle states. Such kinks have been observed in small 
circular dots- "artificial atoms" - because of their special 
circular symmetry More generally, kinks should oc- 
cur in generic quantum dots with no special symmetries 



g. In fact, the data of Ref. Q show evidence for kinks 
in large dots, though they were not the subject of that 
investigation. 

Small circular dots behave much like atoms (hence "ar- 
tificial atoms" ) : the circular symmetry causes degener- 
acy of the orbital levels and so a large spacing between 
allowed energies. In sharp contrast, there is no degener- 
acy in irregular dots: the typical single-particle orbital 
level separation is simply given by A = l/iyV where i' 
is the bulk density of states and V is the volume of the 
dot. Kinks in the evolution of Coulomb blockade peak 
positions may occur whenever the ground state of the 
dot is separated from an excited state with different spin 
by an energy of order A. The interference effects caus- 
ing the separation are unique to each state and change 
upon tuning. In fact, the two states may switch at a cer- 
tain point, the former excited state becoming the ground 
state: such switching corresponds to a kink. Note that 
kinks occur in pairs: kinks in the peaks corresponding to 
the 7V-I-1 transition and that for iV— 1 — > N both occur 
when Egy{N) switches. We see from this argument that 
kinks in the evolution of peak positions with parameter 
is a general property of quantum dots ||] . 

Here we analyze these kinks explicitly in a particular 
limit. Consider a large quantum dot in which the single- 
particle properties are "random" : the statistics of the 
energies follow the classic random matrix ensembles |^ 
and the wave-functions obey Gaussian statistics with a 
correlation function given by the superposition of random 
plane waves ]lO| . The single-particle properties of such 
random systems have been extensively investigated, and 
it has been conjectured, with considerable evidence, that 
these are good models for quantum dots in which the 
classical dynamics is chaotic ||7|,|ll|]. 

To treat the Coulomb blockade, we must consider not 
only single-particle properties but interactions among the 
particles as well. One natural way to proceed is to treat 
these interactions in the basis of self-consistent single- 
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particle states {ipm{^)x<y{s)}^ where m and a are the la- 
bels of orbital and spin states respectively (we neglect the 
weak spin-orbit interaction). In the limit of zero inter- 
action, two electrons occupy each filled orbital state, ex- 
cept for the top level when the total number of electrons 
is odd. Because of the interference produced by confine- 
ment in the dot, the electron density is not smooth but 
rather has small deviations from the classical-liquid re- 
sult. Due to the electron-electron interaction, these devi- 
ations contribute to the ground-state energy, in addition 
to the conventional "classical-liquid" charging energy. If 
the interaction does not change the double-occupancy of 
the levels, one finds that the part of such contribution 
coming from the last doubly-occupied level n is 

^„,„ = y"drdr'[|V„(r)p - {\Mr)\^)]Vscr{r - r') 

x[\Mr')\'~{\Mr')\')]- (1) 

It is appropriate to use the short-ranged screened inter- 
action Vsci- here since the smooth background of the other 
electrons provides screening; (. . .) denotes the standard 
ensemble averaging. If, because of the interactions, one of 
the electrons of that level is promoted to the next orbital 
state, the result (|l|) is modified to become 

itn+i ^ J drdr'[\Mr)\' - (|^„(r)P)]V;cr(r - r') 

x[|V'„+i(r')p-(|^„+i(r')P)] (2) 
± J drdrVn(r)C+i(r)T^scr(r - r')C(r')V'n+i(r') . 

The signs + and — in Eq. (Q) correspond to the singlet 
and triplet states respectively. 

For a large ballistic dot, kpL ^ 1, the lack of correla- 
tion among the random wavefunctions ipn and tpm with 
n ^ m leads to a hierarchy of the matrix elements of 
the interaction [|l2| (here kp is the Fermi wave vector 
of electrons in the dot, and L is the linear size of the 
dot). The first integral in Eq. (||) vanishes in the limit 
kpL oo, and one is left only with the second, exchange 
interaction, term. In the same limit, the exchange inter- 
action term has a non-zero average value and vanishingly 
small mesoscopic fluctuations. The lowest of two energies 
^nn+i corresponds to the triplet state (± — > — ). Neglect- 
ing the small mesoscopic fluctuations of the energies ^„_„ 
and Tt-i-i> O'^s finds the difference between the energies 
of the singlet state formed by doubly-occupied levels and 
of the triplet state with two singly-occupied orbital levels: 

? ^ {U,n) - (C„ + l> (3) 

= 2 fdrdr%,,ir - r') |F(r - r'f = 2V;cr(k-k'). 



Here F{r) = exp(ik'r), with the bar denoting the average 
over the Fermi surface |k'| = kp- 



In the above argument we have implicitly assumed the 
absence of time-reversal symmetry. For i? = 0, the same 
basic argument holds, but the interaction in the Cooper 
channel should be included; this increases ^ making the 
proportionality constant in Eq. (3) larger than 2. 

We thus consider a model jl^ with a single non-zero 
interaction constant ^. This quantity is simply related 
to the usual electron gas parameter (the ratio of the 
Coulomb energy at the mean interparticle distance to 
the Fermi energ y): ^ « (l/A/27r)r,ln(l/r,)A for small 

and in the absence of time-reversal symmetry. As 
increases, the considerations discussed here apply until 
^ exceeds A at which point the Stoner criterion for a 
magnetic instability is approached. For instance, for — 
1 averaging the Thomas-Fermi screened potential over 
the Fermi surface yields ^ = 0.6A in two dimensions. 

The distribution of electrons among the levels depends 
on the single particle level spacing compared to ^. This 
is particularly clear when the total number of electrons 
N is even: the top two electrons can either be in the 
same orbital level at a cost of ^ in interaction energy or 
one can be in level N/2 and the other in N/2 -I- 1 at 
an energy cost of — eN/2- Since the magnitudes 

of both ^ and £n/2+i ~ ^n/2 are of order A, sometimes 
the top level is doubly occupied and sometimes not. In 
the case of double occupation, the state is, of course, a 
singlet; if two sequential levels are occupied, the exchange 
interaction leads to a triplet state. If at most two orbital 
levels are singly occupied, the ground state energy of a 
dot is, then, a sum of three terms: 

£;gr=^ch+ ^n,a+M^ (4) 

(n,(T) filled 

where M is the number of doubly occupied levels. For 
our arguments, the number of electrons is constant and so 
the charging energy Ech is irrelevant. Note that the en- 
ergy (^ is equivalent to the simultaneous filling of two se- 
quences of levels, one of which is rigidly shifted by ^ from 
the other Finally, if several orbital level spacings 

in sequence are small, more complicated configurations 
occur for both even and odd N |l4|jr^ ]. Moreover, the 
problem of the ground-state spin of a mesoscopic system 
becomes very complicated upon approaching the Stoner 
instability p^ . 

As a parameter is varied, the single particle energies 
change and may cause a change in the level occupations 
and so a kink. This is explained and illustrated in Fig. 1. 

The distribution of the kinks in the parameter space 
follows from a random matrix model. We assume that the 
Hamiltonian of the dot follows the Gaussian Orthogonal 
Ensemble in the absence of a magnetic field (GOE, (3 = 
1) and the Gaussian Unitary Ensemble for nonzero B 
(CUE, (3 = 2) 1^. The dependence on the parameter X 
is included by means of a Gaussian process ]T6| , p7t ; we 
consider the process 
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FIG. 1. Mechanism for kinks in the Coulomb blockade peak 
positions: Gaussian process simulation results as a function 
of scaled parameter X (/3 = 2) . (a) 3 typical orbital energy 
levels; note the infrequent avoided crossings, (b) Two sets of 
energy levels, one for singly occupied orbitals (solid) and the 
second for double occupancy (dashed). For a large dot, the 
doubly occupied levels are rigidly shifted from those for single 
occupancy by the interaction ^ (= 0.5 here), see Eq. (^. Con- 
sider an even number of electrons in which the top two parti- 
cles are placed in the levels shown. When the doubly-occupied 
state (first dashed line) has higher energy than the next or- 
bital level (second solid line), the ground state has two par- 
tially occupied levels- a triplet spin state. When these two 
levels cross, the lowest energy in the triplet sector crosses the 
lowest singlet state, (c) Difference in energy upon adding an 
electron, proportional to the Coulomb blockade peak position. 
(Traces offset by 0.3 for clarity.) Crossings of the singlet and 
triplet levels cause kinks in the peak positions (arrows mark 
2 examples). Two adjacent peaks are affected by each change 
in a ground state: they have kinks at the same X and switch 
which orbital level they are tracking. 

H{x) = Hi cos{X) + H2 sm{X) (5) 

where the distributions of Hi and H2 follow the appro- 
priate Gaussian ensemble. Extensive work on parametric 
correlations has shown that the properties of Gaussian 
processes are universal when the energy is measured in 
units of the mean level separation A and the parameter is 
scaled by the typical rate at which energies are perturbed. 



{{de/ dXY)^/"^ JT6| , pT|jT^ ] . For simplicity we restrict our 
attention to kinks caused by a change in configuration of 
the top two electrons when N is even. 

The first quantity to consider is the mean density 
of kinks, pkink- First, because a kink occurs when 
^ = this density is proportional to the 
probabihty of having such a level-separation. Second, 
Pkink must reflect the rate at which the levels change as a 
function of parameter. In fact, it is known that the distri- 
bution of the slopes of the levels, de/dX, is Gaussian and 
independent of the distribution of the levels themselves 
]I6[ . Thus the two contributions are simply multiplied: 

(X for small ^ 

where pi^\s) is the distribution of nearest-neighbor level 
separations, for which the Wigner surmise |^ is an ex- 
cellent approximation, pkink has a strong dependence on 
(3 when ^ is small because of the symmetry dependence 
oi pi^\ and so will be sensitive to a magnetic field. In 
fact, the sensitivity to magnetic field could be used to ex- 
tract experimental values for in quantum dots- a direct 
measure of the strength of interactions. 

Next, an important property is the spacing in X be- 
tween two neighboring kinks. For ^ ^ A, kinks oc- 
cur when two orbital levels come very close and so 
are caused by avoided crossings of single-particle levels. 
Since each avoided crossing produces two kinks, kinks 
will occur in pairs, with small intra-pair and large inter- 
pair separations. The behavior near an avoided cross- 
ing is dominated by just two levels and characterized by 
3 parameters- the mean and difference of the slopes far 
from the crossing and the smallest separation. Wilkinson 
and Austin have derived the joint probability distribution 
of these parameters for Gaussian random processes . 
By expressing the intra-pair separation in terms of these 
parameters and integrating over the joint distribution, 
we find that the distribution of intra-pair separations is 

P{ntra{x) = 2—^ I du exp(-u) (7) 

\ u3/2(l„y^2/^2)-l/2/0p^ P^^- 

For small x, Pintra is linear in the separation x both with 
and without a magnetic field. 

The separation between pairs is usually large for small 
^- avoided crossings with small gap are rare- typically 
many correlation lengths. Hence there is no correlation 
between pairs: Pinter (2;) is Poisson (exponential) for large 
X. Correlation suppresses the probability of two close 
crossings. We make a simple model for this by assuming 
Pinter oc X for X < xq and so approximate Pnter by 
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FIG. 2. Distribution of the separation between near- 
est-neighbor kinks in the Coulomb blockade peak positions. 
The interaction strength is ^ = 0.5. Top: zero magnetic field 
(time-reversal symmetry, l3 — 1). Bottom: non-zero mag- 
netic field (/3 = 2). Insets: tails of the distributions. Both 
the analytic theory for small ^ (dashed line) and the results 
of Gaussian process simulations (solid line) are shown. The 
excellent agreement is remarkable considering that ^ = 0.5 is 
not small and the absence of any fitting parameter. 
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-Pinter (a;) — 

a 



exp(— (a; 
x/xo , 



xo)/a), X > xo 

X < Xq 



(8) 



where C is for normalization. a;o should be of order 1 
in scaled units; we choose it to be the minimum of the 
correlation function of de/dX, xq = 0.85 (0.6) for GOE 
(GUE). The mean density Eq. (||) sets alpha, 
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(9) 



combined with the distribution Pintra- In this way we 
have a parameter free expression for the distribution of 
the separation between adjacent kinks. 

While the above theory is for ^ ^ A, in the exper- 
iments 1^,^ Vg ^ I so that ^ ^ A. Hence we turn to 
numerical calculation to test the range of validity of our 
expressions. Gaussian processes were produced using the 
Hamiltonian (|^) with matrix size 200 over the full inter- 
val X e [0,27r]; the middle third of the spectrum was 
used. Sample energy levels for /3 = 2 are shown in Fig. 1. 

Fig. 2 shows the distribution of kink spacings for 
the experimentally relevant value ^ — 0.5 and f3 — 1, 2. 
Though outside the regime of immediate applicability, 
the theoretical curves agree closely with the numerical 
data. Thus the simple two-level calculation captures the 
main features of the kink distribution for ^ ^ 1 and so is 
adequate for describing experiments in large dots [^. 

In conclusion, through the properties of kinks in the 
Coulomb blockade traces, experiments on quantum dots 



can directly determine the main interaction parameter 
in these dots and so obtain the exchange energy for the 
chaotic electron gas. 
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